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We have investigated the effect of a difference in the tunnelling resistances of the individual nor- 
mal metal-insulator-superconductor (NTS) tunnel junctions in a double junction SINIS device, with 
particular emphasis on the impact on the conductance, cooling and thermometric performance. 
By solving the electrical and thermal equations of the junctions in a self-consistent way, we find 
that asymmetry gives rise to many new features, such as appearance of an excess sub-gap current, 
improved cooling performance, exhibition of negative differential resistance, and improved temper- 
ature range of thermometric sensitivity. Experiments were also carried out to complement some 
of the numerical results. In addition, we studied theoretically and experimentally the effect of a 
finite series resistance, which also causes an excess current in the subgap region, and a suppression 
of the conductance maxima at the gap edge. Experimental results agree well with the theoretical 
predictions. 



INTRODUCTION 

For some time now normal metal-insulator- 
superconductor (NIS) tunnel junctions devices have 
been used as solid state coolers and accurate low tem- 
perature thermometers in the sub-Kelvin temperature 
range [ll-0|. The concept of NIS solid state cooling 
is based on the selective extraction of hot electrons 
from the Fermi tail of the normal metal, using the 
finite energy gap of the superconductor as a filter, 
thereby reducing the average energy of the normal metal 



electrons. Thermometry with NIS ; 



unctions 



, pho nonic 



has 



been used in the study of photonic , _^ , 

electron-phonon [§-12| and electronic 13 1 heat flow, 
and applications in bolometry have been envisioned 
[II, [l3, Ei Cooling has been demostrated for both 
electrons[l,[l6, 17 1 and phonons with high power 

( ^ 50 pW) coolers bein g ab le to cool external payloads 
from 300 mK to 190 mK [li]. In addition, NIS junctions 
have recently been ulitized as accurate charge pumps 
for metrological applications In most of these 

applications, it is beneficial to connect two junctions in 
series to form a SINIS device, essentially doubling the 
performance. 

Traditionally, the performance of SINIS devices has 
been analyzed with two simplifying assumptions - (i) 
the device is symmetric, i.e., the junction resistances 
are equal, and (ii) the normal metal has a vanishingly 
small resistance compared with the tunneling resistances. 
While these assumptions are fair for many devices, one 
should, nevertheless, study the effects of relaxing the 
above assumptions, especially when considering novel 
materials and more complex device processing steps, 
where control over device parameters is more difficult. 

In this paper, we present theoretical and some cor- 
responding experimental results on the impact of asym- 
metric tunnelling resistances, and of a finite normal metal 
series resistance on the conductance, cooling, and ther- 



mometric characteristics of SINIS double tunnel junction 
devices. The evolution of the characteristics as a function 
of increasing asymmetry and series resistance has been 
investigated. A knowledge of these results is not only es- 
sential for a more accurate interpretation of SINIS junc- 
tion based cooling and thermometric experiments, but 
they can also be exploited to optimize the performance 
of such devices. 



SYSTEM UNDER STUDY 

A schematic of a SINIS device is shown in Figure 1, 
where the two junctions are denoted as the left (L) and 
the right (R) junction for convenience, with tunnelling 
resistances Rl and Rfj, respectively. The resistance of 
the normal metal island is Rn- In the case of two probe 
measurements, one also has to consider the effect of the 
probe resistances (Rp), such as the line and the contact 
resistances. These appear in series with the junction re- 
sistances, and can be lumped together with R^ into one 
series resistance Rs = Rn + Rp in the electrical circuit 
equations. The total applied voltage Vr across the device 
is thus 



Vt = Vl + Vr + Vs. 



(1) 



In four probe measurements Vs and Rs are replaced by 
Vn and Rn, respectively. On the other hand, the finite 
resistance of the normal metal island of the device R^ 
means that there is power dissipation in the normal metal 
in the form of Joule heating. This inadvertent heating of 
the device can have an impact on the conductance, cool- 
ing as well as thermometry aspects through the coupling 
of thermal and electrical circuit equations. The Joule 
heating associated with the probe resistances Rp, how- 
ever, does not affect the characteristics of the device, as 
long as all that power is absorbed by the refirgerator far 
from the device. We therefore ignore it in this study. 
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FIG. 1. [Color online] A Schematic of a SINIS device con- 
sisting of two NIS tunnel junction connected in series having 
tunnelling resistances Rl and Rr. Rn is the resistance of the 
normal metal. Ts and Tn are the electron temperatures of 
the superconductor and the normal metal respectively. 



Traditionally, an ideal SINIS junction is analyzed un- 
der the following assumptions - (i) the tunnelling resis- 
tances of the individual junctions are equal (Rl — Rr) 
and (ii) the normal metal resistance is negligibly small 
{Rn ^ Rl,Rr)- This means that the total voltage 
drop across, and the cooling power of an ideal symmet- 
ric SINIS device is simply twice of that for the corre- 
sponding individual NIS junction: Vt — 2Vi = ^Vr and 
Rtot ~ "^Rl = 2Rr. Experimentally, however, these as- 
sumptions are not strictly true as inaccuracies in lithogra- 
phy and variations in local film thickness and morphology 
can give rise to a difference in the tunnelling resistances 
of the two junctions, even when fabricated otherwise un- 
der the same conditions. Moreover, in most actual real 
devices, there is no extra Ohmic contact to the normal 
metal, and therefore a direct determination of the differ- 
ent resistance components is impossible. This means that 
only the total resistance {Rtot — Rl + Rr + Rs) of the SI- 
NIS device is determined form the high bias {Vapp ^ 2A) 
slope of the current / vs applied voltage Vapp curve, and 
the possible effects of asymmetry and series resistances 
are not easily determined. 

For these reasons, it is important to understand how 
the characteristics of a SINIS device evolve as a func- 
tion of the uncertainty associated with the precise deter- 
mination of these individual resistances. Pekola et. al. 
2l| compared theoretically the cooling power of a couple 
of different asymmetric structures at a particular tem- 
perature, and showed that the maximum cooling power 
decreased only a little with increasing asymmetry. How- 
ever, their study was not exhaustive, and a regime where 
cooling can be enhanced can be found, as we show in this 
paper. Also, it is possible that the value oi Rn varies 
independently of the tunnelling resistances, which exper- 
imentally corresponds to changing the material of the 
normal metal island but keeping all other device param- 
eters the same, for example. For the sake of simplicity, 
we consider here the following two distinct cases: 



{i)Effect of asymmetric tunnelling resistance: Rl ^ Rr 
and Rs=0 : If the individual junction tunneling re- 
sistances Rl and Rr differ, the corresponding voltages 
across the individual junctions {Vl and Vr ) are no longer 
identical although the same current fiows through both 
of them. Similarly, the cooling power of the individual 
junctions may vary significantly. 

{ii)Effect of finite series resistance: Rl = Rr and Rs 
y^O : The second goal is to understand the impact of 
the finite series resistances, either in the device itself or 
external. The effects are studied in a symmetric device, 
but could easily be extended to asymmetric cases. 



METHODS 

Theory 

In this work we only consider first-order, single parti- 
cle tunneling theory (22], [23|, and disregard higher order 
two-particle processes |^^24^26] for the clarity of discus- 
sion. The higher order terms become more important 
with transparent barriers (low tunneling resistances) and 
with high resistivity of the normal metal; in many cases 
the junctions are opaque enough and/or the resisitivity 
small enough for the effect to be small. The standard 
result for the current in one voltage biased NIS junction 
is then 



Inis = \ j deNs{e){fs{e) - fN{e + eV)) (2) 

where e is the energy, Ri = ?i/(47re^^|toP-/V'o,AriVo,s) is 
the tunnelling resistance of junction i, determined by the 
area of the junction A, the tunneling probablility j^oP 
and the normal state Fermi-energy densities of states of 
each lead Nq^n/s- fi{^) is the Fermi function in electrode 
i, and Ns{e) is the normalized superconducting quasi- 
particle density of states (DOS). Here, instead of using 
just the equation from BCS theory, we use the expres- 
sion introduced by Dynes et al. [27| to take into account 
non-idealities: 



Ds{e,Ts) 



Re 



iT 



V(e + *r)2-A(Ts)2 



(3) 



where F is the Dynes parameter describing the broad- 
ening of the DOS singularity, and A(Ts) is the super- 
conducting energy gap, with its temperature dependence 
written explicitly. This equation is especially appropriate 
in Al junctions, where it was shown 28| to result from 
environmentally assisted tunneling events (stray photon 
assisted tunneling). One should also point out that, for- 
tuitously, it can also parametrize higher order tunneling 
surprisingly well, including the subgap heating 29] . 
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Since the individual junctions are in series, the same not the same. The expressions for current li and heat 
electrical current flows through both of them, although current Qi through each junction {i = L, R) are thus 
the voltages Vl, Vr differ. However, the heat current is given by the following equations ^J: 



I{VRfL,TN,Ts) = I deNsie,T,)[fsie,Ts) - fNie + eVL,TN)] 

e-Ki J-oo 

= I deNs{e + eVR,Ts)[fN{e,TN)- fs{e + eVR,Ts)], 



(4) 



QLiVLTN,Ts 



e^RL 



de{e + eVL)Ns{e,Ts)[fs{e,Ts) ~ fN{e + Vl,Tn)], 



(5) 



QRiVR,TN,Ts) 



b^Rr 



de{e)Ns{e + cVr, Ts)[/w(e, ^^^) - fs{e + cVr, Ts)]. 



(6) 



The total cooling power of a SINIS device is thus given 



by 



-Ql + Qr, 



(7) 



where the minus sign takes care of the fact that Pt de- 
scribes the total heat power extracted from the normal 
metal. The situation where Pt is positive (negative) thus 
corresponds to cooling (heating) of the normal metal. In 
dynamic equilibrium, this cooling (heating) by the junc- 
tions is balanced by the inflow (outflow) of heat from (to) 
the surroundings, which leads to the equation 



Pt 



th 



+ P[{Pt + I{Vl + Vr)]+I^Rn: 



(8) 

which describes three distinct heating mechanisms: The 
first term on the right gives the direct heating from the 
substrate (phonons at temperature Ttath), which couples 
to the normal metal electrons via electron-phonon inter- 
action [13, H^. B is the coupling strength constant 
which in the present case is given by i? = Sf2 where S 
is the electron-phonon coupling constant and f2 is the 
normal metal electron gas volume. The exponent n = 4- 
6 depends on disorder level, phonon dimensionality and 
film thickness 0, [31-33 1, but is typically n = 5 for good 
metals with thickness > 30 nm. In suspended wires, how- 
ever, phonon transmission out of the nanowire becomes 
the bottleneck for heat flow which leads to B being 
a geometry-dependent constant and n « 2.8. The sec- 
ond term describes the indirect back-flow of heat from 
the superconducting electrodes, associated with, for ex- 
ample, the quasi-particle recombination and backtunnel- 
ing in the superconductor [17|, [3J|. Here (0 < (3 <1) is 
the fraction of total heat deposited in the superconduc- 
tor flowing back to the normal metal. Finally, the third 
term is simply the direct Joule heating of the normal 



r 



metal. Throughout this paper, we take the simplifying 
assumption that the superconductor temperature equals 
the bath temperature, Tg = TBath- This limit is use- 
ful to consider, as it gives the maximum performance 
that a cooler could achieve. A real device can also suffer 
from overheating and/or non-equilibrium effects of the 
superconductor, which degrade the cooler performance 
il[3l[3l. 

In this work, we have carried out numerical calcula- 
tions for two distinct models which we name (1) electri- 
cal and (2) electrical + thermal, for convenience. For the 
first case the calculations are carried out by only consid- 
ering the current and voltage equations ^ and (|4]). In 
this case, all temperatures are equal TBath = Ts = T/v. 
In practice, this limit corresponds to the case where ei- 
ther the normal metal volume is large (large B), or the 
tunneling resistances very high (small Pt and /). In the 
second case, both the electrical [Eqs.([T]) and (U) and the 
thermal equations [ Eq. ([8])] are solved simultaneously 
in a self-consistent way. In the second case, it is then 
possible to determine the value of T/v for a given ap- 
plied bias and TBath , but it is still assumed that Ts = 
TBath- Unless explicitly mentioned, the following typical 
parameter values were used for the simulations : F/A 
=2 X 10-^ /3 =0.03 (from Ref. 0), n=5 and B' =200, 
where we have defined a dimensionless e-p coupling con- 
stant B' = e2A"-2^TB/(1.764/cB)" and Rt = Rl + Rr- 
For typical AlO^; based tunnel junctions with a tunneling 
resistance of 10 kfi, A =220 /j,eV (Al superconductor) 
and E = 1.8 xlO^ (Cu normal metal), the value of B' 
= 200 corresponds to a normal metal electron volume of 
roughly 0.09 fJ,m^^ (for example dimensions = 10 /xmx 300 
nmx30 nm, the normal metal dimensions of our experi- 
mental sample discussed below) . Throughout this paper, 
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FIG. 2. [Color online] Scanning electron micrograph of a SI- 
NIS device with four NIS tunnel junctions. Two SINIS pairs 
of junctions can be connected in such a way that both pairs 
have about the same total tunneling resistance, but one pair 
is asymmetric and the other symmetric. 



for the simulation results, the voltage, current, power, re- 
sistance and temperature are expressed in units of A/e, 
A/ei?T, A^/e^i?T, Rt and Tc, respectively. 



Experimental 

NIS junctions with Al, AlOx and Cu as the super- 
conductor, insulator and the normal metal, respectively, 
were fabricated on oxidized silicon substrates using elec- 
tron beam lithography and shadow evaporation tech- 
niques 0. The thickness of Cu and Al were ~ 30 and 
20 nm, respectively. The insulator was formed by oxi- 
dation of the aluminium surface at room temperature at 
400 mbar for 4 min in pure oxygen. In order to study 
the effect of asymmetry of the tunnelling resistances, 
we fabricated SINIS devices with the geometry depicted 
in Fig. [21 The asymmetric SINIS pair consists of two 
NIS junctions with junction areas 1030x330 nm^ and 
210x330 nm^ (outer junctions), while the corresponding 
symmetric junction pair (inner junctions) have dimen- 
sions 350x330 nm^. The asymmetric SINIS had a total 
tunnelling resistance of 42.4 kl7, quite close to the value 
for the corresponding symmetric SINIS, ~ 44.4 kfi. The 
measurements were carried out using a He^-He'* dilution 
refrigerator with a base temperature of ^ 60 mK. Each 
measurement line had two RC filters, one at 4 K and the 
other at 60 mK, contributing to a total resistance of ~ 
2 kf] per line. In addition, additional microwave filter- 
ing was achieved between the RC filters with the help of 
Thermocoax cables 37|. DC TV curves were measured 



using voltage and current preamplifiers (DL Instruments 
1201 and 1211), and conductance was measured using a 
lock-in measurement with an excitation voltage 8 /iV at 
a frequency 17 Hz. For the measurement of the impact 
of finite resistance, we added external resistors at room 
temperature to the symmetric junction circuit. 



EFFECT OF ASYMMETRIC TUNNELING 
RESISTANCE 

In this section we consider solely the effect of the 
asymmetric tunneling resistances on the SINIS device 
characteristics. We assume Rn =0 such that Rtot — 
Rt — Rl + Rb- First, we neglect the thermal effects 
and try to understand the results in a simple framework. 
Later on, we include the thermal effects. 



Current-voltage and conductance characteristics 

In Fig. [3] (a) we show the evolution of the theoreti- 
cal current - total voltage (/ — Vr) characteristics of a 
SINIS tunnel junction pair with an increasing amount 
of asymmetry in the tunneling resistances, calculated for 
Ts/Tc =0.05. This calculation was carried out in the 
framework of electrical model assuming T/y —Ts, and 
keeping the total tunneling resistance Rt — Rl + Rr 
constant, but varying the values of the individual tun- 
nehng resistances. The case where Rl = Rb = 0.5 Rt is 
the symmetric case, and only values of Rl > 0.5Rt are 
plotted because the results are symmetric with respect 
to the interchange of L and R junctions. It is clearly 
seen that any deviation from the symmetric case gives 
rise to an excess current in the intermediate subgap re- 
gion {bVt = A.. 2 A), whose magnitude increases with 
the asymmetry. Above the subgap region the current is 
reduced with increasing asymmetry. The magnitude of 
the current at eV = 1.5A for Rl/Rt = 0.95 is about 
one order of magnitude higher than in the symmetric 
case, while the increase in current for a smaller asym- 
metry of Rl/Rt — 0.7 is just few tens of percent. The 
differential conductance dl /cIVt curves corresponding to 
Fig. [3{a) are shown in Figure [Sljb) . We see that a large 
enough asymmetry gives rise to an apparent secondary 
quasiparticle DOS "shoulder" between A...2A, with an 
concomitant downward shift of the « 2 A DOS peak and a 
depletion beyond 2A. However, far beyond the gap-edge, 
all the current-voltage and conductance curves merge. 

Let us now try to understand the origin of the excess 
current in the framework of circuit theory. FiglSjc) shows 
the I-V characteristics of a single NIS junction for differ- 
ent values of the tunnelling resistance at Ts/Tc = 0.05. 
Clearly, the shape of the curves are exactly the same, 
but simply scaled by the value of the tunnelling resis- 
tance, as seen from. Eq. Q. By considering a con- 
stant current bias in the region where one of the junc- 
tions is near the gap edge, whereas the other not (for 
example, the dashed line), we can see that the total volt- 
age drop Vt = Vl + Vr is smaller for higher asymmetry, 
or, conversely, for a constant total voltage the current is 
higher, explaining thus the origin of the excess subgap 
current. To understand better why the excess current 
is pronounced around Vr = A.. 2 A, one can study how 
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FIG. 3. [Color online] Calculated (a) current-total voltage and 
(b) differential conductance characteristics of a SINIS device 
at T /Tc ~ 0.05 for different level of asymmetry, calculated 
in the framework of the electrical model. Different plots cor- 
responds to different combinations of Rl and Rr such that 
the total tunneling resistance is constant, (c) Current-voltage 
characteristics of a single NIS junction at T /Tc ~ 0.05 for dif- 
ferent values of the tunnelling resistance. Dashed line gives 
leRr/A = xlO"''. (d) Variation of the voltage of the indi- 
vidual NIS junctions Vl and Vr corresponding to some of the 
plots in (a) as a function of total voltage Vr- 



the individual voltages Vl and Vr vary as a function of 
the total voltage Vr across the SINIS junction, shown in 
FiglSfd). For a low total voltage Vr < A, both junc- 
tions are in the subgap region where the junction resis- 
tances Rji = Vi/I are much higher than the tunnelling 
resistances Ri, but their ratio is the same as the ratio 
of the tunneling resistances Rl/Rt- Therefore, most of 
the voltage drops across the junction with the larger re- 
sistance [Rl in FiglSUd)). With Vt > A, however, the 
resistance of this junction drops dramatically below that 
of the other junction as the gap edge has been reached, 
thereby giving rise to a increase in the current, saturation 
of the junction voltage Vl, and a much faster increase of 
the voltage of the other junction, which has now become 
the high resistance junction (still in subgap). Finally, the 
second junction reaches the gap edge {Vt > 2A), and the 
roles of the junctions are again switched back to the orig- 
inal, with voltage drops increasing as if the junctions are 
simple resistors with zero bias offsets Vl = Vr ~ A. 




FIG. 4. [Color online] (a) Calculated current-voltage charac- 
teristics of a SINIS device at four bath temperatures {T /Tc = 
0.05,0.25, 0.55 and 0.7), and for three different amounts of 
asymmetry Rl/Rt = 0.95, 0.8 and 0.5, calculated without 
taking in account thermal effects. The corresponding con- 
ductance characteristics calculated at bath temperatures of 
Ts/Tc — 0.25 and 0.55 are shown in (b) and (c), respectively, 
(d), (e), (f) are the same, but now including the effect of the 
thermal circuit as well, with B' — 200. 



Next, we consider the evolution of the electrical char- 
acteristics of an asymmetric SINIS device with the bath 
temperature. In Fig. |3] (a) we show the I-Vt charac- 
teristics with three different asymmetry levels at four 
distinct bath temperatures, whereas in Fig. U (b) and 
(c) we plot dl/dVT at two different bath temperatures. 
Clearly, the effects of asymmetry are pronounced in the 
entire temperature regime below Tc- The distinct bump 
in the current-voltage and the conductance characteris- 
tics at Vt = A.. 2 A smears with increasing temperature 
into a broader feature, resembling a shift in V at interme- 
diate temperature range (e.g. at Ts/Tc = 0.25). From 
the conductance curves we observe that the height of the 
conductance maximum (at Vt ~ 2A for low asymmetry) 
decreases and shifts to lower values of Vt as the asym- 
metry increases. The higher the temperature, the larger 
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FIG. 5. [Color online] (a)Measured current-voltage character- 
istics of the SINIS tunnel junction device of Fig. [2] at bath 
temperatures T =0.09,0.2,0.3,0.5,0.7,0.9 and 4.2 K, for the 
symmetric (open circles), and asymmetric (filled circles) junc- 
tion pairs with Rl / Rt ~ 0.875. The lines at 0.7 K (dashed for 
symmetric, solid for asymmetric) are the theory curves calcu- 
lated using parameters corresponding to the sample. The in- 
set shows the results of a measurement at 90 mK in the deep 
subgap region, including the theory curves (solid for asym- 
metric, dashed for symmetric), (b) The measured differential 
conductance characteristics at bath temperatures T — 0.1, 
0.3, and 0.7 K for both the asymmetric pair (filled circles) 
and symmetric pair (open circles). The solid and dashed lines 
are the theoretical curves for the asymmetric, and symmetric 
junction pairs, respectively. 



this shift is: for example at Ts/Tc — 0.55 and for Rl/ Rt 
= 0.95 the relative shift is ^ 40 %. The observed effects 
thus have the implication that if one does not know the 
level of asymmetry, an error in the determination of A 
from the data can follow. Figs. |4Ud)-(f) show the same 
results, but now including the thermal circuit, Eq. [S] 
The general effect of asymmetry is the same, with only 
slight modifications to I-V and dl /dV-V characteristics 
in comparison to the pure electric model, due to self- 
heating and cooling effects. 

In order to verify our computational results, we car- 
ried out experiments to compare the characteristics of an 
asymmetric junction to its symmetric counterpart with 
the sample of Fig. [21 The results are shown in Figure [S] 
(a) for the current-voltage, and (b) for the conductance 



characteristics [38|. In general, the experimental data 
resembles the theoretical plots of Fig. |4] (a) quite well. 
The slight 4 % offset in current between the two curves 
at 4.2 K (where junctions are normal) is simply due to 
the 4 % difference in the total tunneling resistances Rt- 
In agreement with the numerical results, there is clearly 
a voltage dependent excess current in the asymmetric de- 
vice at temperatures below 4.2 K (for example an excess 
of 40 % at 0.3 mV at 500 mK). At the lowest measured 
bath temperature of 90 mK, the predicted hump struc- 
ture in the asymmetric device is clearly visible, shown 
in the inset of Fig. |4] (a). For comparison with theory, 
we have also plotted theoretical curves without free pa- 
rameters at 90 mK (inset) and at 700 mK, with the level 
of asymmetry determined from the SEM image (differ- 
ence in junction areas). The match with experiment is 
very good. The conductance characteristics measured at 
three distinct bath temperatures are also shown in Figl5] 
(b), including the theory curves. The predicted lowering 
and shift of the conductance peak is visible in the exper- 
imental data, especially at the higher temperatures. The 
theoretical curves also match quite well. 



Cooling power characteristics 

The next interesting question is, how are SINIS device 
cooling characteristics affected by the asymmetry? Fig- 
ure inja) shows the calculated bias voltage dependence 
of the total equilibrium cooling power of a SINIS device 
at Ts/Tc =0.25 (where cooling is strongest) for differ- 
ent levels of asymmetry, but keeping Rt constant. We 
see that the cooling power is actually smallest for the 
symmetric junction case, and increases with asymmetry. 
For the highest asymmetry plotted, the enhancement is 
quite significant, a factor of 2.4 at the maximal cooling 
power bias point Vt ^ 1-7 — 1.8A. Fig. [5] (b) shows how 
the cooling power is divided between the two junctions 
vs. Vt in the case where Rl/Rt — 0.8. Interestingly, 
the individual junctions exhibit peak cooling powers at 
different values of the bias voltage, as the individual junc- 
tions reach the gap edge at different times. However, the 
total cooling power is still singly peaked, even though 
there are bias values and bath temperatures, where only 
one junction cools and the other one heats the normal 
metal. 

In order to investigate the cooling aspects further, in 
Fig. IHlc) we show the variation of the maximum cool- 
ing power with bath temperature for various levels of 
asymmetry. These calculations were carried out taking 
into account the temperature dependence of the BCS 
gap, A(T). What we see is that at low temperatures 
{Ts/Tc < 0.4 — 0.5) the peak cooling power increases 
with increasing asymmetry, while at higher temperatures 
the peak cooling power drops with increasing asymmetry. 
The optimal temperature value for cooling [maximum of 
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FIG. 6. [Color online] (a)Variation of the total cooling power 
of the normal metal of asymmetric SINIS junctions as a func- 
tion of applied voltage at Ts /Tc =0.25, with varying asymme- 
try, (b) Voltage dependence of cooling power of the individual 
L and R junctions when Rl/Rt ~ 0.8 (lines with symbols), 
together with the total cooling power (line). Ts /Tc =0.25. 
(c) Variation of the maximum cooling power with tempera- 
ture and asymmetry. The peak shifts to lower temperature 
with increasing asymmetry. The dashed line and the hollow 
symbols correspond to individual L and R junctions when 
Rl/Rt = 0.95. (d) The variation of the maximum cooling 
power (left axis, black symbols) and the temperature of the 
peak in (c) (right axis, red symbols) vs. logarithmic deviation 
from Rl/Rt ~ 1. 
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FIG. 7. [Color online] (a) The variation of temperature of the 
normal metal Tjv as a function of the bias voltage and asym- 
metry at different bath temperatures Ts/Tc =0.05, 0.25, and 
0.55. (b) The variation of the voltage drop across the individ- 
ual high- (L, lines) and low-resistance {R, lines with symbols) 
junctions as a function of the total voltage, calculated at Ts 
=0.05 and taking in account the thermal circuit, for three dif- 
ferent amounts of asymmetry. The markers L and H indicate 
the regions where negative differential resistance is observed, 
(c) Current - voltage characteristics of the individual L (lines) 
and R (lines with symbols) junctions at Ts/Tc =0.05, corre- 
sponding to the the plots in (b). (d) The variation of Tjv as 
a function of the voltage drop across the individual L (lines) 
and R (lines with symbols) junctions at Ts/Tc =0.05. 



Impact of thermal effects: Cooling and negative 
differential resistance 



the peak in Fig. [He)] shifts to lower temperatures, and 
the operational range in temperature (width of the peak) 
narrows with increasing asymmetry. The optimal cooling 
power value keeps increasing up to a very large asym- 
metry Rl/Rt — 0.999, after which the enhancement is 
quickly destroyed, so that in the limit Rl/Rt 1, the 
single junction cooling characteristics are restored [Fig. 
E] (d)]. Fig. [HJc) also shows the variation of the cooling 
power in the individual L and R junctions for the case 
where Rl/Rt = 0.95. It is seen clearly the that junc- 
tion with lower resistance has much higher cooling power, 
but the high resistance junction compensates by heating 
below Ts/Tc = 0.5. 



A more accurate picture of the electrical characteris- 
tics can be obtained by taking into account the thermal 
circuit, Eq. |S1 as well. This now enables us to calcu- 
late the electron temperature of the normal metal T^ as 
a function of the applied bias voltage. Initially, we do 
not include the Joule heating in the normal metal and 
set Rn — 0, later in section □ we also consider the finite 
normal metal resistance. The current-voltage and the 
conductance curves were already shown in Fig. 2] (d)-(f). 
The main influences of the thermal circuit are that the 
excess current feature spreads out to a bit higher volt- 
ages, and the conductance peak shift and reduction is 
not quite as prominent, especially around Tg where cool- 
ing is strong [Fig. 111(e)]. 

For the cooling characteristics, the variation of T^v with 
bias voltage and three different asymmetries are shown 
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in Fig. [3 (a) for Ts/Tc =0.05, 0.25 and 0.55, using the 
scaled e-p parameter B' = 200. At Ts/Tc =0.05, the 
electron temperature first increases with increasing bias 
voltage, because of the finite Dynes parameter heating. 
This heating gives way to the usual strong cooling in the 
vicinity of the gap edge, with T/v reaching its minimum 
at Vt ~ 2A. Curiously, although the cooling power in- 
creases with increasing asymmetry (Fig. |6]), Tjq is not 
much affected. This apparent contradiction appears be- 
cause the cooling power also depends directly on T^v, 
and the asymmetry enhancement effect disappears for 
Tn « Ts- At the intermediate range Ts/Tc =0.25, 
however, the enhanced cooling power is clearly seen as 
temperature reduction at Vr = 1.7A. On the other 
hand, at the highest bath temperatures plotted, asym- 
metry leads to a clear increase of T/v due to the loss of 
cooling power at higher values of Ts, as seen in Fig. [6] 
(c). 

In Figl?] (b) we show the variation of the voltages Vl 
and Vn as a function of total applied bias voltage across 
the SINIS device at Ts/Tc = 0.05, including the ther- 
mal circuit. By comparing Fig. ^h) to Fig. ^d), we 
see that the most distinct feature that appears with the 
thermal effects in high asymmetry devices is the abrupt 
decrease and a subsequent minimum in the voltage across 
the low resistance junction {Rr), above Vr = 2A. This 
region is marked by a box and label H for better visi- 
bility. In order to understand the origin of such a phe- 
nomenon, one has to look at the I-V characteristics of the 
corresponding individual junctions, shown in Fig. [7jc). 
Clearly, the junction with the lower tunneling resistance 
exhibits negative differential resistance (NDR) at 
the high current range, while the corresponding higher 
resistance junction exhibits a much weaker NDR at a 
lower current region around /c-Rt/A = 10"'^. Junc- 
tions with asymmetry Rl/RtO-5- 0.6 (not shown here) 
exhibit no such NDR effect, nor do we observe any NDR 
in the pure electrical model, so that the effect is thermally 
driven. By observing the variation of Tn as a function 
of Vl and Vr (shown in Fig[7| (d)), we notice that NDR 
occurs at regions where there is a sharp change of T/v. 
One should also point out that the observed NDR effect 
is completely stable: the total I-V characteristics never 
exhibit the NDR effect, as one junction always stabilizes 
the other. Also, the NDR effect occurs only in the low 
temperature region, surviving up to Ts/Tc ~ 0.3 for 
asymmetry Rl/Rt — 0.95. A more detailed explanation 
of the NDR effect is given in the Appendix. 

The fact that the strong switching in Tn gives rise to 
the NDR has been investigated further by considering 
the variation of T/v as a function of the applied bias for 
different values of the parameter B' oc BRt, which quan- 
tifies the ratio of heating or cooling by the substrate (e-p 
interaction) to the heat flow through the junctions. In 
practice, this could be changed by changing the normal 
metal volume. In figure[5Ja) we show the variation of T/v 




FIG. 8. [Color online] (a) Variation of Tjv with Vr for 
Rl/Rt =0.95, for (a) various values of parameter B' at 
Ts/Tc =0.05 and /3 = 0.03, and (b) for different values of 
/3=0,0. 03,0. 1,0. 2,0. 3, 0.5,1, with B' = 200. (c) I-V characteris- 
tics of the individual low-resistance junction R corresponding 
to (a). Strong signature of NDR is observed only for small 
values of B' . fi = 0.03. (d) Variation of Tjv with Vr for the 
values of 13 in (b), at Ts/Tc =0.05, B' = 200. 

as a function of voltage drop across the lower resistance 
junction Vr, when Rl/Rt = 0.95, for various values of 
B' . The corresponding current- voltage characteristics of 
the individual NIS junction are shown in[8Kc). Clearly, 
the NDR effect is strong for small values of B' and van- 
ishes when B' > 10^. Similarly, when the backflow of 
heat from the superconductor is increased (parameter /3), 
the dip in the value of T/v decreases (less coohng), result- 
ing in the reduction of the rapid temperature change, and 
the NDR effect, as seen in Figs. [5Jb) and (d). For a dif- 
ferent level of asymmetry, the parameter region in (i?', (3) 
plane where the NDR occurs, changes. 



Thermometry 

For thermometry, a SINIS device is usually current bi- 
ased, so that the thermometric performance can be char- 
acterized by plotting the Vt vs. Tg dependence. This 
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FIG. 9. [Color online] Variation of the total SINIS junction 
voltage, Vt (lines-symbols), for different values of asymme- 
try, as a function of bath temperature Tg, for bias currents 
(a) JeiiT/Alx 10"^, and (c) leRr/^ =1 in the electrical 
model. With thermal effects taken into account, (a) and (c) 
are modified to (b) and (d), where we used B' = 200. The 
open symbols show the individual junction voltages Vl and 
Vr corresponding to Rl/Rt =0.95 (right axis scales). 
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FIG. 10. [Color online] Measured Vt vs. bath temperature 
for both the symmetric (open circels) and asymmetric (filled 
circles) junction pair of the sample, with Rl/Rt ~ 0.875 
for the asymmetric case. The current bias used was 0.7 nA. 
The lines show the theorical curves calculated using known 
parameters of the sample obtained from I-V curves. 



is shown for an asymmetric SINIS device in Figl9l calcu- 
lated using both the purely electrical (Figsl9](a) and(c)) 
and electrical-l-thermal models (FigsE] (b) and(d)), for 
two different values of bias current. Interestingly, for the 
low bias current {IcRt/A = Ix 10""^), the symmetric 
thermometer has the best responsivity, dV/dT. In ad- 
dition, if the thermal circuit (self-heating) is taken into 
account, the thermometer saturates below T/Tc ^ 0.05 
for all values of asymmetry (Fig. [S] (b)). This satura- 
tion is due to the electron-phonon decoupling; T^r does 
not follow Ts anymore, and Vr saturates to a value that 
is determined by the self- heating effects (value of B'). 
Thus, for thermometric applications, it is desirable to 
minimize the effects of the thermal circuit by maximiz- 
ing B' , which can be achieved by increasing the normal 
metal volume, or by increasing Rt- At a bias current high 
enough beyond the subgap regime IcRt/ ^ — 1, shown in 
Figs. Fig. [9] (c) and (d), an asymmetric junction device 
has a better responsivity than the symmetric one at low 
temperatures T /Tc < 0.4, where the symmetric device 
response dies out. However, with the thermal circuit, the 
response again saturates even for the asymmetric junc- 
tions below T /Tc ~ 0.25, leaving only a smaller window 
of improved performance around 0.25 < T/Tq < 0.4 for 
the value of B' = 200 used here. Figures [5] (a)-(d) also 
show the variation of the voltages across the individual 
junctions, showing that the performance degradation at 
low bias currents is caused by the saturation of the low re- 
sistance junction, whereas at high bias the low resistance 
junction gives the performance boost discussed above. In 
conclusion, for high bias currents, a highly asymmetric 
junction can act as a good thermometer over a larger 
temperature range, especially if the thermal circuit ef- 
fects are minimized. This is because one of the junctions 
acts as a good thermometer at low temperatures, while 
the other one has better responsivity at higher tempera- 
tures. 



Finally, we have experimentally checked the effects of 
asymmetry on thermometric performance for the device 
discussed in section „ Fig. [TU] shows the measured Vr vs 
Ts curves for both the asymmetric and symmetric junc- 
tion pairs, for a bias current of 0.7 nA corresponding to 
/ei?T/A ^ 0.13, together with the corresponding theory 
curves of the electrical model. The symmetric theory has 
no fitting parameters, whereas the asymmetric theory has 
only one, the level of asymmetry Rl/Rt ~ 0.875. This 
corresponds well with the ratio of junction areas from the 
SEM image (0.83). It is clear that the general trends of 
the changes due to asymmetry are reproduced, includ- 
ing the improvement in responsivity below T < 0.6 K. 
The discrepancy between theory and experiment at high 
temperatures is not fully understood at the moment. 
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EFFECT OF FINITE SERIES RESISTANCE 

In this section, we consider the effect of finite series 
resistance on the conductance, coohng and thermometry 
aspects of a SINIS device. As already mentioned earher, 
the finite series resistance always models the resistance of 
the normal metal island Rn, but could also include line 
and contact resistances, if the measurement is performed 
in two-probe configuration. In the thermal circuit, only 
Rn appears in the Joule heating term in Eq. [H as the 
Joule heat of the other series resistors is always dissipated 
far from the active device area. For the sake of simplicity, 
we consider here the case of a symmetric SINIS device 
only, Rl = Rb , with results shown for both the electrical 
and electrical + thermal model. Furthermore, for the 
electrical + thermal model, we only include the most 
important resistance Rn, and do not consider the effects 
of other series resistances, which could be easily included 
in the calculations by adding them in the electrical circuit 
equations. 



Impact on conductance characteristics 

While discussing the effect of finite normal metal resis- 
tance on SINIS characteristics we consider two different 
cases: 

(i) In the first case, the value oi Rn is varied keeping 
the total resistance constant (i.e. Rt— Rl + Rr + Rn 
=const.). This situation is typical in actual experiments, 
where the value of Rn remains unknown but Rt is known 
precisely. In Fig[TT] (a) and (c) we show the current- 
voltage and conductance characteristics as a function of 
Rn, calculated at Ts/Tc = 0.05 for the electrical model, 
and in FiglTT] (b) and (d) for the electrical -I- thermal 
model, respectively. It is clearly seen that in both mod- 
els, the subgap current actually increases with increas- 
ing magnitude of Rn, whereas above the gap the cur- 
rent decreases with increasing Rn , until finally all curves 
converge asymptotically at high voltages (not shown). 
This initially counterintuitive behavior can be explained 
as follows: In the subgap region, the resistances of the 
junctions are much larger than values of Rn considered 
here, so that the current is limited by the junctions. In 
that case, increase of Rn corresponds to a decrease of 
Rl+Rr, which leads to increase in the current. However, 
above the gap the junction resistances become compara- 
ble to Rn values considered, and voltage starts to drop 
across Rn- In that case, an increase of Rn leads to a 
decrease of the voltage across the junctions, leading to 
a decrease in current. In terms of differential conduc- 
tance, the main effect in the electrical model (Fig. [TT] 
(c)) is a decrease of the conductance peak value, with a 
concomitant flattening and an outward shift of the peak 
position. If one considers the thermal circuit, the flatten- 
ing of the conductance peak can be arrested as a result 
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FIG. 11. [Color online] Calculated I-V characteristics of a 
symmetric SINIS device vs. Rn, when the total resistance 
Rt~ Rl + Rr + Rn is constant in the (a) electrical and (b) 
electrical -I- thermal case. Ts/Tc = 0.05. The conductance 
characteristic corresponding to (a) and (b) are displayed in 
(c) and (d). 



of the self-cooling and the Joule heating in Rn , but the 
conductance maxima still drops with increasing Rn, 
seen from Fig fTTT d). Thus, similar to the case of asym- 
metric junctions, the net effect of finite Rn is an excess 
current in the subgap and a reduction of current above 
the gap. However, the conductance curves behave in a 
different way both in the subgap and in the conductance 
peak region, as seen by comparing Figs. Sfb) and llll 

(ii) In the second case, Rn assumed to vary indepen- 
dently of i?L + Rr, corresponding experimentally to the 
case where different samples are compared. In FigH^ (a) 
we show the current-voltage characteristics calculated at 
Ts/Tc — 0.05 for the electrical model. Now, the sub- 
gap current is insensitive to the value of Rn, whereas 
the current above the gap is reduced due to the direct 
effect of Rn limiting the current. The addition of the 
thermal circuit gives the same result (not shown), with 
only the shape of the I-V curve changing slightly [as seen 
in Fig. [11] (b).] The behavior of the conductance in 
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FIG. 12. [Color online] (a) Calculated I-V characteristics of 
a symmetric SINIS device vs. Rn, when Rn varies inde- 
pendently of Rl + Rr, for the electrical case. Ts/Tc = 0.05. 
(b) Experimental data (symbols) for conductance curves with 
varying external series resistor (from 4 kQ to 104 kQ) mea- 
sured at r = 0.15 K , together with the theoretical charac- 
teristics (lines) corresponding to the experiment. 



the electrical model is also similar to case (i), with the 
only difference that the high voltage asymptotes natu- 
rally do no converge. This is seen in Fig. [T^ (b), where 
we also plot the experimentally measured characteristics, 
produced by changing an external series resistor. The 
data agrees well with the purely electrical model, clearly 
giving rise to a decrease and shift of the conductance 
maxima at the gap edge. 



Cooling 

The variation of the total equilibrium cooling power 
{Tn = Ts) with total bias voltage Vr = Vl + Vr + Vn 
and series resistance Rn at Ts/Tq = 0.05 is shown in 
Fig. [13] (a) for the case where Rt is kept constant. We 
see that the peak cooling power increases with increasing 
Rn, and shifts to a higher voltage. The power increase 
happens because Rl + Rr has to decrease, thereby in- 
creasing the cooling power, whereas the shift is due to 
the finite voltage drop across -Rat. In the case where Rn 
varies independently, the cooling power naturally doesn't 
change at all, but curves simply shift due to the finite 
voltage drop across Rn (not shown). In addition, if one 
investigates the cooling [Fig. [T3] (b)] and individual NIS 
junction I-V characteristics [Fig. [T3] (c)] as a function of 
Rn in the electrical -|-thermal model with Rt constant, 
one finds again the NDR effect, whose origin is again 
associated with the strong increase of Tn as a function 
of voltage, with the back-bending corresponding to the 
NDR region. The minimum temperature is again not af- 
fected, similar to the asymmetric case [Fig. [7] (a)]. Fig. 
[T31 (d) shows also the single junction I-V curves with Rn 
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FIG. 13. [Color online] (a) Variation of the total cooling power 
of a symmetric SINIS device vs. Vt for different values of 
Rn calculated at Ts/Tc = 0.05 with total resistance Rt = 
Rl + Rr + Rn kept constant, (b) Tn vs individual junction 
voltage for different values of Rn, with Rt constant, (c) I- 
V characteristics of the individual junction with varying Rn, 
for Rt constant, taking in account the thermal circuit, and 
(d) the same but letting Rn vary independently. The points 
also show the result in the purely electrical model for a high 
Rn/Rt = 0.5. 



varied independently, with one I-V curve calculated with 
the pure electrical model and a high value of Rn = Q-5Rt 
also shown. Clearly, no NDR appears in the electrical 
model. The NDR effect is once again fully stable, with 
the stabilization coming from the voltage drop across the 
resistance Rn- 



Thermometry 

The thermometric characteristics, Vr vs. Tg, for the 
first case where the Rt is kept constant, are shown in 
Figs. [M] (a) and (b) for two different bias currents 
leRr/A = 4x10"^ (subgap) and leRr/A = 1 (above 
gap), respectively. Both the purely electrical and elec- 
trical-|-thermal results are shown. As already seen from 
Fig. [11] (a), for a constant / in the subgap region, the 
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FIG. 14. [Color online] Vt vs Ts of a symmetric SINIS device 
for different values of Rn, with a current bias leRr/A — 
4x10""') in (a), and leRr/A = 1 in (b). The total resistance 
Rt = Rl + Rr + Rn is kept constant. The lines show the 
characteristics in the purely electrical model, whereas the lines 
with symbols take into account the thermal circuit. 



voltage drop across the SlNlS pair decreases slightly as 
Rn increases. This is also seen in Fig. [14] (a), where 
it is clear that the effect on responsivity is small. How- 
ever, the situation is reversed for the high bias current. 
Fig. [T3] (b) . Then the responsivity of the thermometer 
is clearly better for higher Rn below Ts/Tq = 0.4, if 
only electrical circuit is considered. On the other hand, 
when Rn varies independently, no performance improve- 
ment or reduction is seen (not shown). At low temper- 
atures the thermometers eventually saturate, with the 
saturation regime depending on the current bias value 
and the thermal circuit (overheating). This overheating 
can negate all the performance gains if B' is not large 
enough, as seen in Fig. [14] (b). 



CONCLUSIONS 



200 % for an extreme asymmetry of Rl/Rt = 0.99. 
If one looks at the maximum cooling power as a func- 
tion of the bath temperature, it shifts to lower tempera- 
tures with increasing asymmetry, and attains an enhance- 
ment of 30 % over the symmetric device. The shift 
of the cooling power maximum to lower temperatures 
is promising for real device operation at lower temper- 
atures. However, the improved cooling power at equi- 
librium does not always translate into a lower minimum 
normal metal temperature. Finally, our theoretical re- 
sults for thermometry show that in the low temperature 
regime T/Tc < 0.2, asymmetry gives typically no bene- 
fits. However, at higher temperatures (using higher bias 
currents) one can take advantage of the different respon- 
sivities of the individual junctions in the asymmetric SI- 
NIS device, to construct a thermometer that is sensitive 
over a broader temperature range. 

We also showed experimentally and theoretically that 
a finite series resistance gives rise to a clear reduction 
of the conductance maximum at the gap edge, and that 
the thermal circuit affects the details significantly. Our 
calculations show that it is also possible to observe the 
negative differential resistance effect in the individual 
junctions in a symmetric device with a high enough nor- 
mal metal resistance. For cooling power, a series resis- 
tance naturally offers no absolute benefits, however, with 
an unknown series resistance, the true cooling power is 
higher than what one would expect for the ideal case. In 
thermometry, there are also no benefits in the absolute 
sense, but an unknown series resistance means that the 
responsivity can be better than what one would expect 
for a current bias above the gap. 
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In this paper, we presented a comprehensive study of 
the effects of (i) asymmetry of the tunneling resistances 
and (ii) additional series resistance on the behavior of 
double junction SINIS tunnel junction devices. Both ex- 
periments and calculations showed that asymmetry re- 
sults in an excess current below the gap voltage and a 
reduction and shift of the conductance peak at the gap 
edge. This effect takes place regardless of the thermal 
circuit characteristics. In addition, our calculations pre- 
dict that it is possible to observe negative differential 
resistance in the individual junctions in a realistic ther- 
mal environment, if asymmetry is high enough. Interest- 
ingly, we also calculate that the total equilibrium cool- 
ing power of an asymmetric device can be larger than 
of the symmetric device with the same total tunneling 
resistance. This enhancement can be very large, over 



APPENDIX 

In order to understand the origin of the observed NDR 
due to asymmetry, let us consider the specific case where 
Rl/Rt = 0.95 and Rr/Rt =0.05. The current- voltage 
characteristics of these individual junctions at two tem- 
peratures Ti and T2 ( where T2 > Ti) are shown in FigfTSl 
Let us consider the situation at a low bias current region 
(shown by marker L). The initial state of the junctions 
L and R at Ti are given II and In respectively. At 
this point we increase the applied bias voltage by a small 
amount, which is sufficient to give rise to an increase of 
Tn from Ti to T2. Since the sub-gap resistance of junc- 
tion L is larger than that of junction R, most of the incre- 
ment in bias voltage drops across the L junction, thereby 
taking it from initial state II to a virtual state VIl- We 
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FIG. 15. [Color online] Schematic representation of the origin 
of NDR. The initial state of the junctions L and R at T\ 
are given II and Ir respectively and the corresponding final 
states by Fl and Fr respectively. When the increase in the 
bias voltage is sufficient to cause a heating of the normal metal 
from an initial temperature T\ to a final temperature T2, the 
junctions goes from the initial state to final state via virtual 
states namely VIl and V2l- The initial, final and the virtual 
states are color coded as green, red and blue respectively. 



identify this state as virtual as we have not considered 
the heating yet. At this point we "turn on" the heating 
which demands that the junction switches to point V2l 
lying on the T2 branch. Once again, we identify V2l as 
a virtual state because if junction L is in V2l, junction 
R has to be in V2ii since they are in series and the same 
current flows through them. However, this is not possi- 
ble, because in that case the sum of the voltage drops 
across the individual junctions exceeds the total applied 
voltage. Thus, the only solution is that for the given 
temperature the total applied voltage is redistributed 
across the individual junctions in such a way that final 
configurations of the L and R are given Fl and Fr re- 
spectively. As seen from the figure, this redistribution 
gives rise to an increase in the current, increase in the 
voltage across i?, a drop in the voltage across L junction, 
and thus leads to negative difi'erential resistance (NDR) 
in the L junction. Using similar arguments, the origin 
of NDR in the R junction in the high current region ( 
corresponding to marker H) can be understood. 
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